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1 Introduction 



In 1876 Smith [16j presented a formula for the determinant of the nxn matrix 
having the greatest common divisor of i and j as its ij element. 
During the 20th century many other results concerning matrices with similar 
structure were published, see for example [71 [T21 [IS]. In 1989 Beslin and 
Ligh |3] introduced the concept of a GCD matrix on a set S, where S = 
{xi,X2, . . . , Xn} C Z"*" with xi < X2 < ■ ■ ■ < Xn and the GCD matrix (S) has 
its ij entry. Since then numerous publications have appeared in 
order to universalize the concept of GCD matrix. For example, Haukkanen [5] 
and Luque [13] consider the determinants of multidimensional generalizations 
of GCD matrices and Hong, Zhou and Zhao [5] study power GCD matrices 
for a unique factorization domain. 

Poset theoretic generalizations of GCD matrices were first introduced by 
Lindstrom [TT] and Wilf [IB]. In these generalizations (P, ^) is a poset, / is 
a function P — > C, S* = {xi, X2, ■ ■ ■ , x„} C P, Xi ^ xj ^ i < j and (5) / is an 
nxn matrix with f{xi A Xj^ 3jS its ij element. These matrices are referred 
to as meet matrices. The papers by Lindstrom and Wilf [IH| arose from 
needs for combinatorics and became possible since Rota [IS] had previously 
developed his famous theory on Mobius functions. Rajarama Bhat [2] and 
Haukkanen [B] were the first to investigate meet matrices systematically, 
presenting many important properties of ordinary GCD matrices in terms 
of meet matrices. In [ID] Korkee and Haukkanen define and study the join 
matrix [S]f of the set S with respect to /, where /(x, V Xj) is the ij element 
of the matrix [S]f. 

During the last ten years the concept of meet matrix has been generalized 
even further in many different ways. Korkee [H] studies the properties of a 
matrix Mg'^'^'^ , which yields both the matrix (S) / and [S]f as its special case. 
A totally different approach is taken by Altinisik, Tuglu and Haukkanen in 
[2], when they define meet and join matrices on two subsets X and y of P. A 
further idea of generalization is presented by Bege [3j as he studies yet another 
GCD related matrix {F{i, {i,j))), where F{m.,n) is an arithmetical function 
of two variables. For present purposes it is convenient to use a slightly 
different notation. For every i G Z"*" we define an arithmetical function /, of 
one variable by 

fi{m) = F{i,m) for all m G Z^. (1-1) 
With this notation Bege's matrix takes the form 

/2((2,1)) 

. /n((n,l)) 

In order to distinguish between this and the numerous other generalizations 



/i((l,2)) ■■■ A((l,n)) 1 
/2((2,2)) ■■■ /2((2,n)) 

U{n,2)) ■■■ /„((n,n)) . 
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of GCD matrices, this matrix is referred to as the row- adjusted GCD matrix 
of the set {1,2,..., n}. This notation also enables us to define row-adjusted 
meet and join matrices. 

Definition 1.1. Let (-P, ^) be a lattice, 5" = {xi, a;2, . . . , a:„} be a finite 
subset of P with Xi ^ Xj ^ i < j and /i,/2,...,/„ be complex- valued 
functions on P. The row-adjusted meet matrix of the set 5" is the n x n 
matrix which has {fi{xi A xj)) as its ij element. Similarly, the 

row-adjusted join matrix ['S']/^,...j„ has {fi{xi V Xj)) as its ij element. 



More explicitly. 



and 



J/lv,/lI 



/i(xiAa;i) /i(a;iA2;2) 
/2(x2Axi) /2(a;2A2;2) 

fn{Xn A Xi) fn{Xn A X2) 



/i(xiVxi) /i(xiV2;2) 
/2(x2Vxi) /2(a;2V2;2) 

/„(x„Vxi) /„(X„VX2) 



fl{Xi A Xn) 
f2{x2 A Xn) 

fni^Xfi A Xn) 



fl{Xi V Xn) 
f2{x2 V Xn) 

fni-^n V Xn) 



;i.3) 



;i-4) 



It turns out that there are some results concerning the matrix (S) fi,...j„ to be 
found in the literature by Lindstrom |11J and Luque [13] . When the notation 
is the same as defined in fll.lj] . these results can easily be apphed to Bege's 
matrix. 

Unlike the ordinary meet and join matrices, the matrices (5')j^^.,,j„ and 
['^]fi,-;U usually not symmetric. There are also many other key properties 
of meet and join matrices that do not hold for row-adjusted meet and join 
matrices. Hence, neither the traditional methods of meet and join matrices 
works in the study of these row-adjusted matrices. 



Remark 1.1. In the case when /i = /2 
{S)f and = [S]f. 



fn = f, we have {S)f^^, 



Remark 1.2. Taking the transpose of a row-adjusted meet or join matrix 
results in a column- adjusted meet or join matrix. Therefore the results con- 
cerning row-adjusted meet and join matrices can easily be translated for 
column- adjusted meet and join matrices using this connection. 

At the end of his paper Bege [5] presents an open problem regarding 
the structure and the determinant of the matrix {F{i, {i, j))). It appears 
that the question about the determinant could be solved using Lindstrom's 
result in [TT]. In this paper we present a more systematic investigation of 
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the structure of ('S')/i,...,/„ and ['S']/^,...j„ in general case. Then by using this 
knowledge we are able to find a different proof for Lindstrom's determinant 
formula and also prove some other results concerning the rank and inverse 
of these matrices. 

2 Preliminaries 

Let (P, :<) be a lattice, S = {xi,X2, ■ ■ ■ , Xn} a finite subset of P and 

/l,/2,...,/n:P^C 

complex-valued functions on P. We also assume that the elements of S are 
distinct and arranged so that 

Xi ^ Xj => i < j. 

The set S is said to be meet closed if x A ?y G for all x,y E S. In other 
words, the structure {S, :<) is a meet semilattice. The concept of join closed 
set is defined dually. 

Let D = {di, (^2, • • • , dm} be another subset of P containing all the ele- 
ments Xi A Xj, i, j = 1,2, . . . ,n, and having its elements arranged so that 

di ^ dj ^ i < j. 

Now for every z = 1, 2, . . . , n we define the function "^nji on D inductively 
as 

^D,/.(4) = /.(4)- E ^A/,(4), (2.1) 

or equivalently 

Mdk) = E "^DjAdv). (2.2) 

Thus we have 

^D,/i(4) = E Md^)tJ-D{dv,dk), (2.3) 

where fin is the Mobius function of the poset (D,^), see jH Section IV. 1] 
and [TTl 3.7.1 Proposition.]. 

Let Eij be the n x m matrix defined as 

" { otherwise! ^^'^^ 

The matrix Ed may be referred to as the incidence matrix of the set D with 
respect to the set S and the partial ordering ^. 

Finally, we need another n x m matrix T = {vij), where 

Vij = {eD)ij'^DjAdj)- (2.5) 

In other words, if S is the n x m matrix having D.fi{dj) as its ij element, 
then T = Ed o S, the Hadamard product of the matrices Ed and S. 
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3 A structure theorem 



In this section we give a factorization of the matrix ('5')/j_...j„, which then 
enables us to derive formulas for the rank, the determinant and the inverse 
of the matrix (5*)/^^.. ^,^. 

Theorem 3.1. We have 

{S)f,_f^=TEl = {EnoE)El. (3.1) 

Proof. By ({22]), (I23D and (j^SD the ij element of (5)/,,...j„ is 

m 

f{xiAXj)= J2 ^A/.K) = EMife*A/.(eD),fc, (3.2) 

dv^Xi/\Xj k=l 

which is the ij element of the matrix TEj). □ 

Remark 3.1. It is possible to define row-adjusted meet and join matrices 
/!,...,/„ and [X,Y]f-^^^ j^ on two sets X and Y by [{X,Y)f-^^ j^)ij — 
fi{xiAyj) and ([X, y]/i,...j„)ij = fi{xi\/yj). It would be possible to generalize 
Theorem 13.11 for these matrices, but the methods used in the proofs of the 
other theorems do not work in this general case. 

Remark 3.2. In the case when the set S is meet closed Theorem 13.11 also 
provides an effective way to calculate all the necessary values ^ sjii^j) as 
follows. In this case D = S and both Eg and T are square matrices of size 
n X n. Since Es is also invertible, from equation fl3.ll) we obtain 

T = (5)/,...j„(i?J)-\ (3.3) 

which gives the values of "^sjA^j)- Here the matrix E^ is the matrix as- 
sociated with the zeta function (s of the set S (see pfl p. 139]), and thus 
the matrix {Eg)~^ is the matrix of the Mobius function of the set S and has 
^s{xi,Xj) as its ij element. 

The following example gives a solution for the first part of Bege's problem. 

Example 3.1. The row-adjusted GCD matrix of the set S = {1,2, ... ,n} 
is the product of the matrices T = {vij) and Eg, where 

^""'^'^ = otherwise ^^'^^ 



and 



n 



esh'^sjXJ) = ies)ij M^)f^ T = (eshifi * /^)(i), (3.5) 



k\j 

where * is the Dirichlet convolution and fi is the number-theoretic Mobius 
function. It should be noted that here the notation F{i,k) = fi{k) is not 
only convenient but also enables the use of the Dirichlet convolution. 
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4 Rank estimations 



In this section we derive bounds for rank (5')/^^..,j„ in the case when the set 
S is meet closed. The rank of meet and join matrices or even GCD and LCM 
matrices has not been studied earher in the hterature. 

Theorem 4.1. Let S be a meet closed set and let k be the number of indices 
i with £,j^{xi) = 0. Then the following properties hold. 

1. rank {S) f^,...j„ = iff fi{xi A Xj) = for all i,j = l,...,n. 

2. If k = 0, then rank (S')j^^...j^ = n. 

3. Ifk> 0, then 

n-A;<rank(5)/,,...j„ <n-l. (4.1) 
Proof. 1. Follows trivially. 

2. By Theorem 13.11 we have 

rank (5)^,...,;„ = rank (TEJ) . (4.2) 

Since in this case the matrices T and Es are both triangular square 
matrices with full rank, the claim follows immediately. 

3. Since multiplying with the invertible matrix Eg does not change the 
rank, we have 

rank (S')/j....j^ = rank T. (4-3) 

To obtain the latter inequality we only need to note that since at least 
one of the diagonal elements of T equals zero, the rows of T cannot be 
linearly independent and thereby T cannot have a full rank. On the 
other hand, the n — k rows with nonzero diagonal elements constitute 
a linearly independent set, from which we obtain the first inequality. 

□ 

In the case when the set S is meet closed and fi = ■ ■ ■ = fn = f (that is 
in the case of ordinary meet matrix) the question of the rank becomes trivial. 
Namely, the matrix [S) f can be written as 

{S)f = EsAE^, (4.4) 

where A = diag (\l/5j(xi), '^s,f{x2), ■ ■ ■ , \E'5j(x„)), see [2] Theorem 3.1.]. Now 
by the same argument as in the proof of Theorem 14.11 we have 

rank {S) j = rank A = n — k. (4-5) 

The following two examples show that the bounds in Theorem 14.11 are the 
best possible under these assumptions. They also show that a large value of 
k may indicate a large decline of the rank of the row- adjusted meet matrix, 
but not necessarily. 
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Example 4.1. Let xi = Xj A Xj for all i, j = 1, . . . ,n, which implies that xi 
is the smallest element of 5* and the set iS'\{xi} is an antichain. Now the set 
S is clearly meet closed, and for every z = 2, . . . , n we have 

<fsJAx^)=fiix,)-fi{x^). (4.6) 

If i > 1 and we set fi{xi) = fi{xi), then the ith column of T becomes the 
zero vector and thus for every i > 1 we may reduce the rank of the matrix 
{S)fi,...j„ by one. Therefore if the first diagonal element of T is not zero, 
then rank (S) = n — k. 

Example 4.2. Let (P,^) = A/s and = F as shown in Figure [H Let 

f2ix2) = hixi) = fsixs) = fiixa) = Uixi) = h{xi) = f^{x^) = 1 (4.7) 
and fi{xj) = otherwise. Simple calculations show that '^sj2{x2) = 1 7^ 0, 
^sjAx,) = ^s,dxs) = ^s,.uix^) = "^sj.ix^) = 0, (4.8) 
and thereby k = A. But on the other hand we have 



iS)fl,-,fn - 





















1 








1 


1 


1 


1 


1 


1 








1 


1 


1 











1 


1 



and clearly rank (S) fi,...j„ = 4. 




Figure 1: The lattice A/s and the choices of the elements of the set S. 
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5 Determinant formula 



In this section we present a determinant formula for the matrix (<S')/j_...j„ 
when the set S is meet cfosed. This theorem is almost the same as that 
presented by Lindstrom [TT]. It is possible to use the Cauchy-Binet equality 
to obtain a determinant formula for in general case. Since it is 

similar to the case of usual meet matrix, we do not present it here. 

Theorem 5.1 ([H], Theorem). If the set S is meet closed, then 

n n 

det(5)/,,...j„ = n^5,/,(2:i) = n fi{xj)ns{xj,Xi). (5.1) 

i=l i=l Xj<Xi 

Proof. Since the set S is meet closed, we have D = S. Then the matrix 
Es is a lower triangular square matrix having every main diagonal ele- 
ment equal to 1. The matrix T is a lower triangular square matrix with 
^sjil^i)) ^5,/2(^2), • • • 5 ^5,/„(a;n) as diagonal elements. Thus det Es = 1 
and by Theorem 13.11 we have 

n 

det(5)/„...j„ = detT = n^5,/.(xO. (5.2) 

2=1 

The second equality follows from 02.31) . 

□ 

Remark 5.1. The original theorem by Lindstrom [TT] is slightly more general 
since it does not require the assumption Xi ^ xj ^ i < j . As he states, the 
rows and columns of ('S')/i,...j„ can always be permuted in a way that does 
not change the determinant but makes the matrix [S) f^,...j^ to fulfill this 
condition. 

The following example gives a solution to the second part of Bege's prob- 
lem. 

Example 5.1. For the row-adjusted GCD matrix on the set S = {1, 2, . . . , n} 

we have 

n n f i\ 

det({l, 2, . . . , n})^,.„,;„ = n ^sjXi) = 11 E [-] =X{{h* 

i=l i=l j\i \] / 1=1 

(5.3) 

6 Inverse formula 

In this section we study the inverse of the matrix {S) fi,...j„ when the set S is 
meet closed. A formula for (5")^^^ in general case could be obtained with 
the aid of meet matrices on two sets and the Cauchy-Binet equation. We do 
not, however, present the details here. 
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Theorem 6.1. // the set S is meet closed, then the matrix {S)f^„,,j^ is 
invertible iff sjii^i) 7^ for alii = l,...,n. Furthermore, in this case 
the inverse of is the n x n matrix B = {hij) with 



hj = ^fis{xi,Xk)Okj, 

k=j 

where the numbers 9jj, . . . , 9nj are defined recursively as 



(6.1) 



1 



7kj 



k-1 



ifk = j. 



7— tE ^ku"^ s,h{xu)Ouj ifk> J. 



(6.2) 



Proof. The first part follows directly from Theorem 15. 1[ To prove the second 
part we use Theorem 13.11 and we obtain 



(6.3) 



In order to obtain the ij element of the matrix (15)^^^ we only have to 
ascertain the ith row of {Eg)~^ and the jth column of T~^. As stated in 
Remark 13. 2^ the matrix {Eg)~^ is the matrix associated with the Mobius 
function of the set S. Therefore its zth row is 



... fIs{Xi,Xi) fis{Xi,Xi+i) ... fIs{Xi,Xn) 



(6.4) 



Now let G = {6ij) denote the inverse of T. By multiplying the jth row of 
T with the jth column of G, we obtain 



(6.5) 



Further, the multiplication of the kth row of T and the jth column of G 
results in 



^ku"^ S,h{Xu)Ouj = 0. 



(6.6) 



Thus we obtain 06.21) . and (16. ip follows when we multiply the matrices G and 

□ 



7 Formulas for row-adjusted join matrices 

In this section the results presented in previous sections are translated for 
row-adjusted join matrices. The proofs of these dual theorems are omitted 
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for the sake of brevity. Row-adjusted join matrices (or even row-adjusted 
LCM matrices) have not previously been studied in the hterature. As stated 
in Remark II. 2^ the study of column-adjusted join matrices can easily be 
reverted to the study of row-adjusted join matrices via taking the transpose. 

Let D' = {d[, d'2, ■ ■ ■ , d'^,} be a subset of P containing all the elements 
Xi y Xj, i, j = 1,2, . . . ,n, and having its elements arranged so that 

d[^d'^^i< J. 

For every i = 1,2, . . . ,n we define the function "^'^i j_ on D' inductively as 



or equivalently 



Thus we have 



/,(4)= Y: (7.2) 



^^'./.(4)= E MOt^D'idid'^), (7.3) 

where fiD' is the Mobius function of the poset {D', ^), see [T71 3.7.2 Propo- 
sition.] . 

Let E'^, be the n x m' matrix defined as 

" { otherwise. ^'^'^^ 
Finally, let T' = {v[j) be the n x m' matrix, where 

^:, = (e^')..^D',/.K)- (7.5) 

Theorem 7.1. 

[S]f.,...,u = '^'(K'f- (7.6) 

Theorem 7.2. Let S be a join closed set and let k be the number of indices 
i with vl'^, j.{xi) = 0. Then the following properties hold. 

1. rank [S*] /,,...,/„ = iff fi{xi V Xj) = for all i,j = l,...,n. 

2. If k = 0, then rank ['S']/^,...j„ = n. 

3. Ifk> 0, then 



n — k < rank [5']/^,...j„ < n — 1. (7.7) 
Theorem 7.3. // the set S is join closed, then 

n n 

det[5]/i,...,/„ = n*5,/,(^») = n E fiixj)f^s{xi,Xj). (7.8) 

i=l i=l Xi^Xj 
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Theorem 7.4. // the set S is join closed, then the matrix [S]f-^^,,,j^ is in- 
vertible iff' ^'5. j^(xj) 7^ for all i = 1, . . . ,n. Furthermore, in this case the 
inverse of is the n x n matrix B' = (b'^j) with 

b',j = j2l^s{xk,Xi)e',^, (7.9) 

k=l 

where the numbers 9jj, Oj-ij, • • • , O'lj are defined recursively as 



^kj = { __ . ^ . , _ . (7-10) 
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